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Abstract 

We examine the occurrence of Bose-Einstein condensation in both nonre- 

lativistic and relativistic systems with no self-interactions in a general setting. 

A simple condition for the occurrence of Bose-Einstein condensation is given. 

We show that condensation can occur only if q > 3, where q is the dimension 

associated with the continuous part of the eigenvalue spectrum of the Hamil- 

tonian for nonrelativistic systems or the spatial part of the Klein-Gordon 

operator for relativistic systems. Furthermore we show that the criterion for 

the appearance of the Meissner-Ochsenfeld effect is closely connected with 

that for the appearance of Bose-Einstein condensation. 
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I. INTRODUCTION 



One of the most interesting properties of a system of bosons is that under certain con- 
ditions it is possible to have a phase transition at a critical value of the temperature in 
which a macroscopic fraction of the bosons can condense into the ground state. This was 
first predicted over 70 years ago for the ideal nonrelativistic Bose gas ||T],3 and is nowadays 
well known to happen if the spatial dimension D > 3. (See || for the case D = 3 and Q 
for general D.) The analogous phenomenon for ideal relativistic gases has also been stud- 
ied and the criterion on the spatial dimension has been found to be the same as that for 
nonrelativistic gases 

Of course it is of great interest to see what influence different types of interactions have 
on the condensation. In nonrelativistic theories in certain models it is known what happens 
if interactions are included ||. In relativistic theories a detailed study of A</> 4 theory at finite 
temperature and density has been given in flat |j,nj as well as in curved spacetime [pjj . 

Rather than examine self-interacting fields, or the interactions among different quantum 
fields, a simpler problem is to study what happens for a quantum field under external con- 
ditions, where by external conditions we mean for example gravitational or electromagnetic 
background fields or simply boundary conditions imposed on the field. Several different 
situations have already been considered. Let us only mention the analysis in the presence of 
boundaries and static gravitational fields |T2|-[20|| and in the presence of (mainly constant) 
mae netic fields 0-0 ■ 

One of the main questions in all these kinds of considerations is, whether or not Bose- 
Einstein condensation can occur. In the above literature rather involved analysis has been 
done in order to arrive at the conclusion whether or not Bose-Einstein condensation can 
occur. Recently, the authors gave a very simple criterion to decide whether a system might 
condense or not JHJ. Here, we plan to give more details of the calculations and to apply the 
criterion to a large class of examples. This will include the examples treated in the above 
mentioned literature as well as some new cases. Furthermore we will show how the criterion 
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might be used to decide if a system reveals the Meissner-Ochsenfeld effect or not. 

It is worth emphasizing precisely what is meant by BEC since different definitions are 
possible. For the free Bose gas in three or more spatial dimensions the specific heat has a 
non-smooth behaviour at a critical temperature which signals a phase transition. In terms 
of the quantum field theory the phase transition may be interpreted as symmetry breaking 
with the scalar field developing a non-zero expectation value at the critical temperature. 
Associated with this phase transition is a sudden growth in the occupancy of the ground 
state. In the present paper we will adopt this as our definition of BEC. It is possible to 
have a build-up of particles in the ground state without a phase transition occurring. In 
this case there is no unique way to identify a temperature associated with the build-up of 
particles in the ground state. For the particular case of a constant magnetic field in three 
spatial dimensions it was shown by Rojas |32| that although there was no phase transition 



there could still be a build-up of particles in the ground state which could be interpreted as 
BEC. A similar situation occurs for bosons confined by a harmonic oscillator potential [|34| . 



II. CRITERION FOR THE OCCURRENCE OF BOSE-EINSTEIN 

CONDENSATION 

A. Relativistic quantum field theory 

Let us consider a complex relativistic scalar field which may interact with background 
electromagnetic or gravitational fields but which is otherwise free. We restrict our attention 
to an ultrastatic spacetime of the form M. = R x E with metric 

ds 2 = dt 2 — Qij{x)dx l dx^ (2.1) 

and with the field obeying any boundary conditions in the spatial directions. The action 
functional for the complex field $ will be chosen to be 



S = J dt J da x {(D»<$>y(D^)-m 2 &<$>-U [ 



x 



(2.2) 



da x is the invariant volume element on the manifold S. In order to have a scalar field 
exclusively under external conditions, the functions Uq(x) and U\{x) may depend on the 
background gravitational or electromagnetic fields, but are independent of the scalar field. 
They are also assumed to be independent of the time. The kinetic term is the usual gauge- 
covariant derivative 



(2.3) 



In an ultrastatic spacetime, finite temperature and density are easily incorporated in the 
Euclidean time formalism [j35|-|3"7|1 . It has been shown, that the partition function Z may be 
represented in the form (see for example ||38|| ) 



Z = J [d$][d$t]exp{-£} 



(2.4) 



where 



S = J J da? 

o s 



$ f + ie(A -ifi)& 



x 



$ - ie(A - 



+ |D$| 



+ (m 2 + Ui(x))&$ + U (x) + J$ + J+$t 



(2.5) 



In order to decide whether or not Bose-Einstein condensation will occur we need the charge 
in the excited states. It may be obtained directly from the effective action. Eliminating the 
dependence on the sources J and Jt by 

SW , 



sw , 



5 J] IftJ' 

with the background fields $, the effective action is defined by 

r\n,$,&] = w\ti,j,ji] 



(2.6) 



- J dr J da x [J$ + J f $ f ] 



(2.7) 



o s 



with 



W[n,J,f] = -lnZ[/i, J, J*]. 



(2.8) 



It is known, that a minimization of the effective action is equivalent to a minimization of 



the Helmholtz free energy ||25|| . In terms of the effective action the charge is then defined by 

1 dT , 



Q 



(2.9) 



We may now proceed with the analysis of the effective action. It is straightforward to 
show that 



r[/i,$,$t] = + -lndet(/ 2 ^. 
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(2.10) 



where we used DeWitt's condensed notation [ft9|| . I is an arbitrary unit of length introduced 
to keep the argument of the logarithm in fl2.10| ) dimensionless. Specializing to the case of a 
background static magnetic field only, and choosing the gauge 



A = 0, V • A = 0, 



(2.11) 



we have 



1 



lndet(/ 2 ^) = r+ + r_ 



(2.12) 



where 



r± = ^lndet i 1- 







dr 

+m 2 + U 1 (x)]} 



Te/i -D 2 



(2.13) 



For the calculation of (|2.13| ) we will use the zeta function scheme ^y^T[. In this scheme 
one defines 



r ± = -^(o) + ^C ± (o)iu/ 2 



(2.14) 



with the generalized zeta functions 



oo 



(2.15) 



j=— oo N 



Here, X^ N are the eigenvalues of the fluctuation operator (see fl2.13| )) given by 




(2.16) 



with the eigenvalues <jn of the spatial part of the Klein-Gordon operator, 



(-D 2 + U x {x) + m 2 )f N {x) = a N f N (x), 



(2.17) 



with a complete orthonormal set of functions /at(x). It is easily seen that T + = T_, so that 
we skip the index ± in the notation and consider ( + . 

Now all the definitions and preparation to explain our criterion for the appearance of 
Bose-Einstein condensation have been given. For the moment we do not want to consider 
any specific situation but want to assume a quite general structure of the eigenvalues Xjn 
or, more specifically, of the energy E N = a^. As we will see, to decide if Bose-Einstein 
condensation appears, the only relevant information on En is about its continuous part. 
Let us assume that splits into the sum of a discrete part (p is just a set of labels 
for the discrete part of the spectrum), and a continuous part which we can deal with by 
imposing box normalization. The box will be taken to have sides Li, L q , for some q and 
thus we write 




(2.18) 



In the limit Li — > oo the starting point for the analysis of £(s) is 



oo 




(2.19) 



Doing a resummation in I, this is equivalent to 
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x j dt r-^e-^-^ ■<>'< 

o 

(4vr)^ r( S ) U H 2 



EE 



( 2 



- 



8+1 



xK&L^tfly/a*) (e^ 1 + e'W) , (2.20) 
where we introduced the zeta-function of the discrete part of the spatial section, 

a=oo = E«)~ s - ( 2 - 21 ) 



The condition for Bose-Einstein condensation to occur is that fi must reach a critical value 
He set by the lowest eigenvalue in the spectrum, 



,, 2 — rr d — 



(2.22) 



If the charge Q in the excited states, Eq. (|2.9|) , remains bounded as \i — > fic, then Bose- 
Einstein condensation occurs, because for the total charge large enough, it is not possible 
to accommodate it all in the excited states. If Q is not bounded as fi — > fic, then any 
amount of the total charge can reside in the excited states and Bose-Einstein condensation 
will not occur. We therefore need to look at the behaviour of (d/dfi)({0) and (d/dfi)C{®) 
as fi — > fie- 

Using the definition (|2.9|) and representation (|2.20 ) for ((s), one immediately finds 



(4tt) 



i=i p 



g+l 
2 



2 fir 



( 



x a — e 



) 



[2.23) 



The convergence of the sums is defined through the behaviour of the MacDonald functions 



for large arguments [45 



K v {z) ~ ^e- z {1 + 0{z- 1 )} . (2.24) 

It is clear that in the limit /i — > /i<7 all but the contributions coming from p = are finite. 
Thus in the leading approximation, neglecting finite pieces, we write 

rd\ q/2 °° 

^Tr^V^v 7 ^). (2.25) 

We have allowed for the possible degeneracy of the ground state by introducing the degen- 
eracy factor do- It is seen clearly that in the limit fi — > fie — V^o the charge remains finite 
for q > 3; thus for q > 3 Bose-Einstein condensation occurs, and for q < 2 Bose-Einstein 
condensation does not occur. 

The detailed behaviour for \i — ► \iq is most clearly extracted using the Mellin-Barnes 
integral representation of the exponential in Q2.25| ), 

c+ioo 

e~ v = — [ da T(a)ir a , (2.26) 
2m J 

c—ioo 

with $lv > and c 6 1, c> 0. Using this in Eq. (|2.25| ), one finds 

/ m\ 2 c+ioc 



x 



^ - //r a r Q C« (a + |) , (2.27) 



where, in order to allow for interchanging the sum and the integral one has to impose 
a > 1 — q/2. (r(s) is the Riemann ^-function, which is analytic in s except at s — 1 where 
it has a simple pole with residue 1. Closing the contour to the left of the rightmost pole 
then gives the following leading behaviour : 



q = 0, pole of order one at a = 1, 



- fc) = (2-28) 



? = 1, pole of order one at a = 1/2, 



S 



dpVi ( He 



) 



1/2 



(2.29) 



q = 2, pole of order two at a = 0, 



^0^2 

2vr/? 



(2.30) 



As mentioned, for g > 3 no divergent contribution results. The restriction q > 3 includes a 
large number of previously known results, often established by long and detailed calculations, 
as special cases. These and some new examples are summarized in Sec. 3. 



We will consider a nonrelativistic field theory described by the complex Schrodinger field 
$ whose action functional is 



da x is the invariant volume element on the D-dimensional Riemannian manifold E. We will 
consider S to be compact. (In the case E = M D we will impose box normalization with 
the infinite box limit taken.) Ui(x) is an arbitrary time independent potential. As before, 
D = V — ieA is the gauge covariant derivative, with A a vector potential describing a 
background magnetic field and |D$| 2 is short for g t3 DjQ where gr^ is the Riemannian 
metric on E. 

The theory described by (|2.31|) has the local gauge invariance 



B. Nonrelativistic quantum field theory 




(2.31) 



A -> A + W , 



(2.32) 



which gives rise to a conserved current, and a conserved charge which is given by 
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Q 



daJ&\ 2 . 



(2.33) 



We will deal with the charge rather than the particle number to mirror the relativistic case 
considered above as closely as possible. The conserved charge is dealt with by introducing a 
Lagrange multiplier fi, which is the chemical potential. After a rotation to imaginary time 
r, the partition function is again expressed in the form ( |2.4| ), with 



+-^|D$| 2 + f/ 1 (x)|$| 2 
2m 

-e/i|$| 2 + J$ + JT^l . 



(2.34) 

The field $ has been coupled to a complex source J as in (|2.5|) . The effective action can 
now be defined as in Sec. 2.1. In place of (|2.10 ) we find 



1^,$,$*] = +lndetZ 



d 



^2.35) 



where $ is the background Schrodinger field. The second term in (|2.35 ) has arisen from 
performing the functional integral in (|2.4|) over the Schrodinger field. 



In order to regularize the determinant which appears in (|2.35| ) we will again use the (- 
function method. This time we will let /W(x) denote the eigenfunctions of D 2 + t/i(x) : 



/jv(x) = a N f N (^) 



(2.36) 



(These /W(x) differ only by a trivial rescaling from the /at(x) used in Sec. 2.1.) The /at(x) 
are seen to be stationary state solutions to the Schrodinger equation for whatever boundary 
conditions are imposed. The eigenvalues are the energy levels for the first quantized 
system. The set {f N (x)} is assumed to be complete and orthonormal. 

Because the functional integral in ( |2.4| ) extends over all fields periodic in the imaginary 
time coordinate r with period j3 = T -1 , the eigenvalues of the operator appearing in fl2.35p 



arc 
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XjN = 2m jT + cttv — e/i 



(2.37) 



where j = 0, ±1, ±2, .... The generalized (^-function is given by (see [2.15|) 



j=—oo N 



and we define 



In det I 







1 



— -ey,-— D z + [/i(x 



<9r 



2m 



(2.38) 



-C'(0) + C(0)lnZ . (2.39) 



Using the summation formula in [42], it is easy to show that 



C(s) = Y^(cr N -efi)- s + ( T (s) 



(2.40) 



where we have defined 



Ct(s) 



N 



EE 



^ e -n/3(o-jv-e/Lt) 



v > N n=l 



(2A1) 



The first term in ( p. 40 ) has no explicit temperature dependence, and corresponds to the 
zero-point energy contribution to the effective action. We will ignore it in what follows. It 
disappears if normal ordering of the Hamiltonian is adopted. 

The term in Cr(s), which we may call the thermal ^-function, is easily shown to have the 
following values : 



Ct(0) = , 

&(0) = - i 1 - e-^-^] . 



(2.42) 
(2.43) 



N 



These last two results show that the second term in (|2.35|) has a simple relation to the 
thermodynamic potential Q defined by 



Q = T^ln [1 -e 



•/3(o-jv-e/i) 



(2.44) 



N 
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(We could have adopted Q as our starting point; however, our criterion is most simply ex- 
pressed using the generalized C-function, and in addition we wished to parallel the relativistic 
calculation.) 

We will again assume that ctn splits up into the sum of a discrete part <7p, and a contin- 
uous part which is dealt with by box normalization as in ( |2.18| ) with E 2 N replaced by a at. In 
the large box limit, the thermal ^-function becomes 

y T q/2-s ~ e -n(3(a*- ef i) 

Ct(s) = (4tt)*/2 t(s) \ ^ nW- ' (2 ' 45) 

As in the relativistic case, the lowest mode <Tg plays the crucial role in determining whether 
or not Bose-Einstein condensation can occur. We therefore split the sum over p in (|2.45 ) 



by separating off the lowest mode for special treatment. The critical value of the chemical 
potential is given by 

efi C = a = 4, (2.46) 

and we wish to study the behaviour of the charge as \x — > fic- If we write 

(t(s) = C¥\s) + C { / 0) (s), (2.47) 

where 

1 } (4vr)^ 2 Y{s) ^ n^/ 2 -° 1 } 



is the lowest mode contribution (with degeneracy do) , and QT (s) is given by (|2.45 ) with the 



sum over p restricted to non-zero values, it is easy to see that because the argument of the 

d 

exponential in ( |2.45| ) will always be negative, even for fi = fi c , Ct^'(O) anc ^ ~q~ Ct^'(O) 
remain finite. Bose-Einstein condensation is therefore determined from a knowledge of 

cfw. 



From ( gjD we have C {0) (°) = and 



W (4vr)9/2 J 2^ n i +q /2 ■ 

n=l 
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The charge is given by 



where terms which remain finite as \x — > \ic have been dropped. Using ( 2.49|) we have 

Q^^.edod- E ng/2 ■ ( 2 - 51 ) 

\ ' n=l 

At /i = the sum appearing in this expression for the total charge diverges for q < 2. This 
gives us the necessary and sufficient condition for Bose-Einstein condensation to occur. For 
q > 3, Bose-Einstein condensation does not occur. 

We can obtain more detailed information on how Q diverges as \x — > fic as before by 
making use of ( 2.26Q . It is easily seen that (|2.51|) becomes 

\ q/2 pc+ioo 



Qifi ^ ^ = i^r (£) Loo daT{a) 



x ( — ) ( R (a + q/2) . (2.52) 



This result may be used to show that 



Q( M ^ MC )^^L ( g = ); (2.53) 
He - H 

Qin^^-^edoKTinc-H)- 112 (9=1); (2.54) 

Q(^c)^^ln(— ) (* = 2). (2.55) 

Only the leading part of Q which diverges as \x — > \xc has been shown in these expres- 
sions. The way in which Q diverges as /x — > He is seen to be the same for relativistic and 
nonrelativistic systems. 

III. APPLICATION OF THE CRITERION TO EXAMPLES 

In this section we will consider several examples of an ideal Bose gas under possible exter- 
nal conditions and we will see that the conclusion whether or not Bose-Einstein condensation 
might occur can be drawn very easily using our criterion of section 2. 
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Let us start with the free ideal Bose gas in a (D + l)-dimensional Minkowski space-time. 
Then the eigenvalues of the Klein-Gordon operator ( |2.17| ) are 

El = k 2 + m 2 , keR D . (3.1) 

In the non-relativistic case we have 

E k = k e m D . (3.2) 

In either case we see that q = D, and conclude that BEC can only occur for D > 3. 
This agrees with conclusions of Refs. for the relativistic gas, and Refs. [^3|,[|3| for the 

non-relativistic gas. 

There are many possibilities of restricting the Minkowski space by imposing boundary 
conditions in one or more directions. For example imposing Dirichlet boundary conditions 
in one direction, fx(xi = 0) = /jv(#i = L), the energy eigenvalues eq. ( |2.1^ ) are 

E 2 nrk = k 2 +(^y + m 2 , keR D -\ «GN, (3.3) 



for the relativistic field, and 



1 f , /7T7l\2' 



E ^ = 2m~ r + Vl) ) ' k G R ' n G N ' (3 ' 4) 
for the non-relativistic field. In this case we find q = D — 1 and therefore BEC is only 
expected for D > 4 (corresponding to q > 3). If we impose Dirichlet boundary conditions in 
p of the spatial dimensions then q = D —p and BEC would only be expected for D > 3 + p. 
The choice of Dirichlet boundary conditions is not important here, and periodic or Neumann 
boundary conditions would lead to the same result. The restriction D > 3 +p holds for any 
choice of boundary conditions which results in a discrete spectrum in p spatial dimensions. 

Another important case where the energy spectrum contains a discrete part is found 
when there is a constant external magnetic field present. Starting with a single component 
constant magnetic field B one encounters the well known Landau levels giving 

E 2 n % = (2n + l)eB + k 2 ± + m 2 , ne N , k± E R D ~ 2 , (3.5) 
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for the relativistic field, and 



E 



n,k i 



2m 



(2n + l)eB + k] 



(3.6) 



in the nonrelativistic case. In both cases the eigenvalues are degenerate with degeneracy 
(SeB)/{2ir) with the (formal) volume S of 1R 2 . The single component magnetic field results 
in q = D — 2, and condensation can occur only for D > 5. This was shown for the relativistic 
field in [E7LEH] and for the non- relativistic field in ||24j| . In particular BEC is absent for D = 3 



22 



In the general case, a magnetic field in D spatial dimensions is characterized by 5 inde- 
pendent components, where D = 26 or 25 + 1 depending upon whether D is even or odd 
f26| , |38|1 . It is straightforward to show that 



for the relativistic field, and 



(3.7) 



3=1 



E 



n,k i 



2m I ^ y 

\ 3=1 



rij + l)eBj + k 



2 

-L C ' 



for the non-relativistic field. In both cases n.,- G and k± G 



nD-21 



(3.8) 



. The degeneracy is 



Y[ P j=i e BjL 2 j-iL 2 j/ (2tt) where Bj j — 1, . . . ,p represent the independent field strengths and 
p < 5. (Lj is the length of the box in the j th direction with Lj — > oo assumed.) In either case, 
q = D — 2p, so that every component of the magnetic field reduces the effective dimension 
by 2. The criterion for BEC is D > 3 + 2p in agreement with the explicit calculations in 
Refs. 

Another example is to consider a field on the space S p x ]R D_P where S p is the p- 
dimensional sphere. In this case we have 



E' 2 - = n(n + p - 1)gT 2 + k 2 + £a~ 2 + m 2 



n,k 



(3.9) 



for the relativistic field, and 



E t = 

n,K 2m 



n(n + p- l)a~ 2 + £a~ 2 + k 2 
15 



(3.10) 



for the non-relativistic field, a represents the radius of the sphere, and £ is a dimensionless 
coupling constant which is related to a possible coupling of the scalar field to the scalar 
curvature. We have q = D — p here, and BEC requires D > 3 + p. In the case of greatest 
physical interest we have p = D = 3 which corresponds to the Einstein static universe. BEC 
does not occur in this case as shown recently in Ref. [J0[ . 



In the examples we have discussed so far the energy spectrum has been explicitly known 
so that a direct evaluation of the thermodynamics has been possible. Our next example is 
to consider a case where the spectrum is not explicitly known; namely a single component 
constant magnetic field in a cylindrical box of radius R, where the field satisfies Dirichlet 



boundary conditions f38[. The eigenvalues for this problem are given in implicit form only 
and the equation reads 

'n 1 1 



Ml 'j- ^n + l;^eBR 2 j =0 (3.11) 

with n E N, 7 = (2/eB)[E N - eB(n + 1) - kj_), k\ e R D ~ 2 , and M(a,b;c) the confluent 
hypergeometric function |j45| . The zeroes of ( |3.11| ) are a discrete set of values for 7; thus 
q = D — 2 and condensation occurs for D > 5, a criterion derived without any problem. 
The exact analysis of the problem would be impossible due to the lack of an explicit form 
for the zeros of the confluent hypergeometric function. 

A similar situation to the last example occurs if we consider a gas confined by an arbi- 
trarily shaped cavity of dimension p in a space of dimension D. It is impossible to write 
down any explicit form for the energy eigenvalues. If the confinement is achieved by impos- 
ing Dirichlet boundary conditions say, then the portion of the energy spectrum coming from 
the cavity would be discrete and q = D — p. A more complete analysis of this problem will 
be given elsewhere [fl6] . 



As our final example we consider a non-relativistic gas in a harmonic oscillator potential. 
If we assume that 

f/(x) = im^^ 2 , (3.12) 
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with p < D, then 



E 



n,k 




(3.13) 



are the energy levels. We have q = D — p here so that D > 3 + p is required for BEC. In 
the case of greatest physical interest D = 3 it can be seen that BEC does not occur. At this 
stage it is worth repeating the comment in the last paragraph of Sec. 1. We are interpreting 
BEC in the sense of symmetry breaking with an associated phase transition, such as that 
which occurs for the free Bose gas with D — 3. This does not occur unless D > 3 + p. We 
are not saying that there cannot be a build-up of particles in the ground state unassociated 
with a phase transition as the temperature is lowered. 

This concludes the list of examples for the relativistic and non-relativistic Bose gas. As 
we have seen, known examples and in addition new examples can be dealt with very easily. 

IV. APPEARANCE OF THE MEISSNER-OCHSENFELD EFFECT FOR A 
GENERAL HOMOGENEOUS MAGNETIC FIELD 



Let us consider the situation where a general homogeneous magnetic field is present. The 
energy eigenvalues then read 



with the limit Lj = L — > oo taken and Vj = eBj. p is the number of non-zero components 
of the magnetic field as discussed in the last section. The degeneracy for each Landau level 



A. Relativistic scalar field in an external magnetic field 




(4.1) 



is (L 2 /(27r)) p rij=i v j- (This degeneracy corresponds to do introduced in Sec. 2.) Thus for 



the zeta function we have 




(4.2) 



oo 
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For the magnetization 



we need (d/dBi)((s). (The magnetization and the choice of units is described in the Ap- 
pendix.) Using the same kind of steps as for the derivation of (|2.20|) , we find 



s 



Vj=l 



Z 2 7T 2 1 j = 1 J=1 ni) ___ jnj)=0 



e(2n; + 1) 



2 / V°p 



+(C - -I*), (4.4) 



where we defined 



QM= Yl (2^ + l)K)- s , (4.5) 
m,...,n p =0 

and abbreviated to (/i — * — //) the terms obtained from replacing // in the second set of terms 
in (|4.4|) with — In order to decide whether or not the Meissner-Ochsenfeld effect takes 
place we can consider the limit Bj — > and see if a nonvanishing magnetization results. This 
is equivalent to the generalization of the Meissner-Ochsenfeld effect from D = 3 dimensions 
PD| to arbitrary dimension D. When Bi — > this removes one of the discrete degenerate 
Landau levels and replaces it with two additional continuous labels. (See eq. (4.1).) The 
effect of this is to increase q by 2 to q e ff = q + 2. Our general analysis showed that the 



gas could only condense for q > 3. Thus if q e ff > 3, or q = 1, 2 there may be interesting 
behaviour when a magnetic field is present even though strict BEC cannot occur. The case 
q — 1 in the limit Bi —>■ can be compared with the gas for Bi = with q = 3; the case 
g = 2 in the limit -B,; — > can be compared with the gas for Bi = with q = 4. In either 
case we wish to study the behaviour of the magnetization as fi becomes close to the critical 
value fi c . As argued in our criterion for BEC, the leading term comes from the smallest 
energy eigenvalue and it reads 

a-q/2 ( p \ 



^ r 9/2 e -/9J(Mo-M) j (4. 6 ) 



oo 

X 
1=1 



which is very similar to the behaviour of the charge, Eq. (|2.25 ). We see that for q > 3 the 
sum is convergent, even for fi = fj,Q, and thus M^Bi — > 0) = for that case. However, for 
q — 1,2, we find using once more (|2.26|) . 



q=l 



MiiB, - 0) = ( fl ^ ) (4.7) 



M,(5i - 0) = --^-^ ]ln(//c-//). (4.8) 

2 2 7T2 y i=1 y 

For g = in the limit 5j — >• the effective dimension gets increased to q = 2 and thus no 
condensation can occur. Thus /i < /ic, the sums are convergent and though Mj(_Bj — > 0) = 
in that case in the given approximation. 

The magnetization laws, Eqs. ( |4.7| ) and ([Of), strongly depend on how fi — > for 



— > 0. It is clear, that the charge in the continuum above the lowest Landau level will 
condense into the ground state in the limit B^ — > 0. Let us use the notation Q gr (fi — > /ic) 
for that charge. This is however exactly the charge given in Eqs. ( |2.29p and ( ^.3U| ), and 
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the limit Mi(Bi — > 0)/Q gr (n — > He) is seen to be finite. If we introduce the charge density 
Pgr = Q gr /L D , it reads for q — 1, 2, 

M ^ - g) = (4 .9) 

Using the universal knowledge of the charge Q gr in the ground state as a function of the 
critical temperature [I2,I3|, 



D-r 



Q gr = Q M - ^- 
with the total charge Q, and total charge density p = Q/L D , we arrive at 



(4.10) 



Mi{Bi -> 0) 



2^ /? ( 1 -(I 



(4.11) 



with /ic = wm 2 + X^=ij as the magnetization law for q — 1,2. For g = 1 in three 



spatial dimensions (D = 3) the result reduces to the one of ref. p0 |, and we agree with their 
conclusions on the Meissner-Ochsenfeld effect. 



B. Non-relativistic scalar field 



The energy levels for a ^component constant magnetic field are 



En ^± (2nj+1) ^ + ± £ (4 , 2) 

j=l j=2p+l V 3 7 

where rij — 0, 1, . . ., and Lj = L — > oo as before. The degeneracy is L 2p Yl^=i e ^j/(^ n )- The 
effective action is determined as described in Sec. 2.2, and the magnetization found as in 
(4.3). We can again argue that q = 1,2 are the interesting values as Bi — > since in this 
limit we obtain q e ff = q + 2 = 3, 4 for which BEC is possible. The lowest state contribution 
is of greatest interest in this limit and we find 

e / m \ q ^ 2 f eB 
Mi(Bi -> 0) ~ -— 



2m \2np J \ 2vr 

n=l 
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as the leading contribution. (Note that // c = ^ 2~^j=i e -% here.) Apart from a numerical 
prefactor, this is the same as in the relativistic case in (4.6). For q > 3 the sum converges 
and we find Mj — > as B[ —>■ even if fi = fi c . The sum occurring in ( (4.13j ) is the same as 
that occurring in the ground state charge (2.51).. If we write p gr = Q(p ~ /i c )/L D as the 
ground state charge density we find 

Mi(Bi ^ 0) ~ . (4.14) 



This is essentially the same as the result found by Schafroth |22[ for D = 3. 



V. CONCLUSIONS 

In this article we considered the general setting of a quantum field, relativistic or nonrel- 
ativistic, without self-interactions under general external conditions. We tried to answer the 
question whether or not the system might undergo a Bose-Einstein condensation, where by 
Bose-Einstein condensation we mean a phase transition associated with a build-up of parti- 
cles in the ground state. Our main result is that the crucial feature governing Bose-Einstein 
condensation is the dimension q associated with the continuous part of the eigenvalue spec- 
trum of the Hamiltonian for nonrelativistic systems or the spatial part of the Klein-Gordon 
operator for relativistic systems. In either case Bose-Einstein condensation can only occur if 
q > 3. Having this criterion at hand, many of previous results were obtained easily in section 
3. In addition some new applications were given. The criterion was shown by studying the 
lowest eigenvalue of a differential operator. This relates to the idea of the effective infrared 



dimensions and finite size effects studied earlier |^7[ . 

Furthermore, we applied our ideas to the appearance of the Meissner-Ochsenfeld effect 
for a general homogeneous magnetic field. With the help of our criterion a discussion of this 



effect is greatly simplified and special cases are easily recovered [p2| |30| , |3~l| 1 . 

The extension of our method to study Bose-Einstein condensation in theories with self- 
interaction 0,|IUJ is of obvious interest. Recently we have shown how this problem can be 
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tackled very efficiently using ^-function methods [11]. It is very likely that the method 
described in the present paper can be extended to interacting field theory. In particular it 
is possible to define an effective field theory describing the lowest modes as in ref. |47| . 
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APPENDIX: A NOTE ABOUT UNITS 



As in classical electromagnetism, the magnetization M causes an effective current density 



which changes the effective field strength. This was noted by Schafroth |22| in the non- 
relativistic case, who argued that because the radius of the classical orbit for the charged 
particles was much greater than the average interparticle separation (whose scale is set 
by the particle density), the effective field strength should be identified with the average 
microscopic field. In the opposite case, where the interparticle separation is much greater 
than the classical charged particle orbit, then the correct procedure is to treat the particles 
as point dipoles. 

An important point is that Schafroth used cgs units for his electromagnetic field con- 
ventions. In this case (see Jackson ^8[) the magnetization gives rise to an effective cur- 
rent density Jm — cV x M which must be added to the current in the Maxwell equation 

47r Air 
V x B = — J. The resulting Maxwell equation may be written as V x H = — J where 
c ' c 

H = B-4ttM. (Al) 

An important point is that the factor of 4tt which enters here has nothing to do with the 
fact that we are in three spatial dimensions. It is simply a consequence of the fact that cgs 
units have been chosen. In SI units we make the changes 



B -> J— B 

V Mo 
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and (|A1|) is changed to 



H= — B-M . 

Ho 

In this case there is no factor of 4n in front of M. 



(A2) 



Our reason for making such a fuss over the units is that Ref. p7p8| has altered to 
read 

Dtt d / 2 



H = B 



-M 



(A3) 



r(i + D/2) 

in D spatial dimensions. This expression first appears in the work of May who gener- 
alized the original calculation of Schafroth |22| from 3 to D dimensions. The -D-dependent 



factor in (|A3|) was noted as the area of the (D — l)-dimensional sphere by May. However, 
it should be clear that the factor appearing in front of the magnetization M is due to the 
choice of units, nothing more. (For the purist, we note that all of the above can be rewritten 
in terms of differential forms with the exterior derivative d in place of Vx.) 

There is another simple way of deducing the effective magnetic field which bears out our 
interpretation. The energy stored in the magnetic field is 



W = — / dcr x B ■ B 



in cgs units If B is constant, then we have simply 

1 



(A4) 



W 



-VB 1 



57T 



(A5) 



with B = |B|. This term should really be added to the free energy F computed from the 
partition function or thermodynamic potential to form 



F+—VB 2 . 

8tt 



(A6) 
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The magnetization is identified by the usual thermodynamic expression E9] to be 



and therefore we have 



dF T V 

' — (B-AirM). (A8) 



dB 8tt 

In our calculation, M = |M| has arisen from one-loop quantum effects. In the absence of 



quantum effects we would have simply the term in ( |A8|) which involves B. This suggests 



that we define the effective field strength as in (|A1|) . If SI units are adopted the factors of 

— in ( |A4| -[A"6|) are replaced with , and a repeat of this simple argument identifies ( |A2|) 

8n — — 2/i 

as the effective field strength. 

In relativistic quantum field theory it is conventional to use neither cgs nor SI units, but 
instead Heaviside-Lorentz rationalized units. In these units, where the lagrangian density 
for electromagnetism is —-F^F^, ( |A6| ) becomes 

F T = F+ l -VB\ (A9) 

and we identify 

H = B — M (A10) 

as the effective field strength. (The difference between Heaviside-Lorentz rationalized units 
and cgs units is that B cgs = V^tt-Bhl, with a similar scaling relation for the vector gauge 



fields A^. See Schweber [|5(J for example.) We adopt Heaviside-Lorentz rationalized units 
in our work, as usual in quantum field theory. 

We will now show how the effective action formalism may be used to calculate the 
magnetization for a general applied magnetic field in a space of arbitrary dimension. The 
full effective action may be written as 

T = S cm + S + f (All) 
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where 



■%,„ = > / <io, ^F l3 F l '-j: xt A i } (A 12) 



is the electromagnetic field action in Heaviside-Lorentz rationalized units. is the magnetic 
field tensor in a Riemannian space £ with the time taken to be imaginary with periodicity 
(3. We assume that the magnetic field is static, but otherwise arbitrary at this stage. Jg Xt 
is the externally applied current responsible for setting up the magnetic field. S represents 
the contribution from the background scalar field, or Schrodinger field. T represents the 
contribution from quantum effects of the matter fields. 

We can settle the question of the correct expression for the magnetization in D spatial 
dimensions by working out the Maxwell equations for the magnetic field with quantum effects 
included. We simply need to evaluate ST when Ai or F^ is varied, and set the variation 
equal to zero, since this will give the effective field equation. Variation of (|A12[) leads to 



SS em = pjda x {VjF ij ~ JL} SAi . (A13) 

Here Vj is the covariant derivative computed using the metric on E. We will define the 
ground state contribution to the current density by 

SS=-f3 da x J l gmmid SAi , (A 14) 

7s 

as in Ref. [fj"2]| . We now need to compute ST. In the actual calculation of T, the result is 
expressed in terms of Fij rather than directly in terms of A\. However, it is easy to prove 
that 

St (M5) 



5 A J \ SFij 



If we define the current induced by quantum effects J md by |42 



J '-^~hm ] ' (A16) 



ST 

then the effective Maxwell equation obtained from = is 

SAi 
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VjF lJ — J* xt + J l ound + J l md . (A17) 



It is possible to rewrite (|A17|) in a form which resembles the result in the case of D = 3. 



This is done by defining the tensor Hy by 

jjij = pij + 2T- — . (A18) 
SFij 

As we will show in a moment, this is the analogue of the vector H in the usual D = 3 case. 
With the definition in ( |A18| ), we can rewrite ( |A17 ) as 



V^' = JL + j; ound • (A19) 

In order to see that our definitions correspond to the usual ones for D = 3, described in 

the earlier part of this appendix, restrict attention to D = 3. The magnetic field vector B 

with components B l is defined in terms of the field strength tensor Fy by Fy = €ijkB k , or 

equivalently by B % = -e^ k Fjk, where e^j. is the Levi-Civitta tensor. In the same way we 

can define the vector H with components H l in terms of the tensor Hij by Hij = eijkH k or 
1 



H % = -e^ Hjk- Contraction of both sides of ( |A18|) with e^, and using 



results in 



5T 1 . .j, ST , . 

w„ = ? m • (A20) 



Hb = Bt+T iw- (A21) 



The thermodynamic potential Q is given by 

f = pn . (A22) 
If we define the thermodynamic potential density to by 

fl = [ da x oj , (A23) 

then we have 
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Hi = Bi — Mi , (A24) 

where 

^ = ~ (A25) 
defines the components of the magnetization. 
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